Jaw.wa functions in Passarino-Veltman reduction 



D. Bardin, L. Kalinovskaya, V. Kolesnikov and W. von Schlippe 

December 19, 2009 

Abstract 

In this paper we continue to study a special class of Passarino-Veltman functions J 
arising at the reduction of infrared divergent box diagrams. We describe a procedure of 
separation of two types of singularities, infrared and mass singularities, which are ab- 
sorbed in simple Co functions. The infrared divergences of Cq's can be regularized then 
by any method: photon mass, dimensionally or by the width of an unstable particle. 
Functions J, in turn, are represented as certain linear combinations of the standard Dq 
and Co Passarino-Veltman functions. The former are free of both types of singularities 
and are expressed as explicit and compact linear combinations of logarithms and dilog- 
arithm functions. We present extensive comparisons of numerical results with those 
obtained with the aid of the LoopTools package. 
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1 Introduction 



In the standard Passarino-Veltman reduction [L\ of 4-point box functions with an internal 
photon hne connecting two external lines on the mass shell there appears an infrared and 
mass singular Dq function (see, for example, p]). A typical example of these diagrams arising 
in the calculation of one-loop EW corrections to // —>■ ZZ{ZA) processes was considered 
in [3], where a universal approach to the calculation of such diagrams was proposed. 

In this paper we describe how this approach works for t — > bfif[ and fif[ — > tb (/i is a 
massless fermion) Charged Current (CC) processes. For these processes one meets eight such 
box functions, four direct and four cross ones. Cross boxes are trivially derived from direct 
ones by a permutation of arguments. Boxes for i decays are related to those of t decays, 
see Ref. So, it is sufficient to consider only one pair of boxes shown in Fig. [H 
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Figure 1: Remaining (left) and J^^y^ (right) functions. 
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Using the standard Passarino-Veltman reduction it is possible to derive relations between 
infrared and mass singular functions DQ{—ml,—m'f,—m'f^,—m'^,Q^,T'^-,0,rnb,Mn',md) and 
Co{—m'^, —■ml,T'^]rnd,0,mi,) and an infrared finite but mass-singular auxiliary function 
•^AwiQ'^y T^; rub, rrit, rrid, mu, M^) and another Co{—m^, — m^, Q^; M^r, nid, 0) with mass singu- 
larity. These basic relations, exact in masses, are: 



JAwiQ iT ■,mb,mt,md,mu,Mn') 
{Ml + Q')Do{ 



-m 



6) 



-m 



t 5 



m^,Q ,T ■,0,mb,M^,md) 



+Co{-ml, -ml, md, 0) - Co(-m^, -ml, T^; m^, 0, mb), 



J, 



WAiQ iT ■,mt,mb,mu,md,Mn'] 



{Ml + Q^)D,{ 



-m 



by 



-m 



t 5 



m'l, Q'\ T^; M^, mt, 0, m„) 



+Co{-ml, -ml, Q"^; 0, m„, M^.) - Co{-ml, -ml, T^; mt, 0, m„). (1) 



Let us emphasize that we have changed the ordering of mass arguments of J{^^ as compared 
to Dq. For they are ordered into two pairs of heavy (6, t) and light {d,u) quarks such 
that the first mass in each pair corresponds to the fermion coupled to the photon, leading to 
the appearance of a contribution logarithmically singular in this mass. 
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The great advantage of basic relations ([T]) is the following. The complex object Dq, con- 
taining an infrared divergence, is excluded in favor of an explicitly computed J function and 
simple objects, Co functions, whose infrared divergences can be regularized by any method: 
by a photon mass, by dimensional regularization or by the width of an unstable particle. 
Examples of the latter Cq functions, regularized by the width, may be found in Ref. [1]. 

We use the standard SANC definitions: 



The paper is organized as follows. 

Section 2 is devoted to the calculation of one of these functions — J'^^^ = J. From here 
on we omit indices of J since the list of arguments uniquely defines its type. 

In Section 3 we present a similar calculation of the J function for the process ud — tb. 
The direct AW and WA functions are defined by: 

JlwiQ^^ T'^i '^d, fUu, rrib, rrit) = 

(Q^ + Ml,)Do{-ml, -ml, -m^, -ml, T^; 0, ma, M^, mb) 
+Co{-mf, -ml, M^, m^, 0) - Co(-m^, -m^, T^; m^, 0, m^), 

JwAiQ"^^ ^u, rrid, mt, mb) = 

(Q^ + Ml,)Doi-ml, -ml, -ml, -ml, , T^; M^, m„, 0, mt) 
+Coi-ml, -ml, 0, mt, M^) - Co(-m^, -ml, T^; 0, mt). (3) 

Again, we limit ourselves to the presentation of function J'^^r■ 

In Section 4 we briefly discuss the J functions for the t channel process bu — > td. 

For all processes we take the limit of vanishing light quark masses. The mass of the 
quark which is not coupled to the photon may be set equal to zero, while that for the quark 
coupled with the photon develops a mass singular logarithm. We keep logarithmic terms and 
neglect quark masses everywhere else. This approximation results in different expressions for 
J functions for the three channels under consideration, and this is why the derivation must 
be presented for three channels separately. 

Every Section ends by a numerical comparison with results obtained with the aid of the 
LoopTools package [5] for zero width and IR regularization by infinitesimal photon mass. 

Section 5 contains a short introduction to the FORTRAN packages, which realize the 
calculation of "doubly subtracted" J functions (see Sections 2-4 for their definition). 

In Section 6 we present our conclusions. 

2 Calculation of the J function for t bud decay 

2.1 Representation in the form of a triple integral 

The basic definition of the function J reads: 



Q' = {Pl+P2)', 



U'={P2+P,)'. 



(2) 



iiT^J{Q^, T^; mb, mt, md, m„, M^) = fi' 



4-n 



/ 




(4) 
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where 

do = q"^ , di = (q + Piy + ml, 

d2 = {q + Pi+P2? + Ml, ds = {q-p^f + ml. (5) 

By standard Feynman parametrization introducing variables x, y, z, as is shown in Fig. [1], 
one can pass to a (3 + n)-tuple integral over x, y, z and over the internal momentum q. In 
n-dimensional space we have: 

In our case we have = Lz, p = zk^y, a = n = 4, therefore Eq.([6]) becomes: 

In terms of Feynman variables the denominator takes the form: 

D = (io(l — z) + dizyil — x) + d2zxy + d^z{l — y). (8) 
From expression ([8]) we derive: 

D = (f — 2zqkxy + Lz — ie, (9) 

where the variable L and the vector k^y are: 

L = (Q^ + Ml-)xy+{pl + ml){l-x)y+{pl + ml){l-y), 
kxy = PA{l-y) -Piyil- x) - {pi+p2)xy. (10) 

Since —ie is an infinitesimal addition, it is possible to replace it by —iez and redefine L 
and D: 



L = L — ie, 

D = q^ - 2zqk.xy + Lz. (11) 

The triple integral over Feynman parameters may be expressed by the same Eqs. (13), (16)- 
(17) as given in Ref. [3]: 

11 1 



J= I dx I ydyN^y / dz _ , (12) 
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where we have neglected the hght quark mass which does not lead to a mass singularity, 
and changed the notation of the other masses as follows: 



rrib 
rrit 

rrid 



nil , 
m4 . 



The ingredients entering Eq. (1121) are 



xy 



and 



L 

p2 

i2 



y 



-2k^yP2 = -2mlxy + Ny , 
T\l-y) + yml + Q^ + ml, 



-mlx^y'^ + Nyxy - Ty , 
P^xy — ie, 

T'^yiX -y) + fnjy + ml{l - y). 



2.2 Integration with respect to z 

The integration with respect to z is straightforward: 



dz 



(L - Zklyf 



{^xyY 



ln(L - kly 



-lnL+ 



L — kly_ 



(13) 



(14) 



(15) 



(16) 



2.3 Integration with respect to x 

In Eqs. ( !T2|) . (fT6ll we perform, first of all, a change of variables: 

xy = X , ydx = dx , 

hence the ingredients become: 



The key identity is 



xy 

k'^ 
'^xy 

L 



—2m\x + Ny , 
—m^x^ + NyX — Ty . 
P^x — ie . 



N = -e ■ 



(17) 



(19) 



it allows the integration by parts in full analogy with section 2.3 of Ref. [3]: 

y y 



I{y) = j dUV = UV - j dVU, (20) 

^5 

where we have introduced an infinitesimal parameter 6, because both parts in Eq. (1201) diverge 
separately. 
Let 

L* = L-kl^. (21) 

be a quadratic trinomial in x: 

L* = Alx^ + Blx + Cl = Al{x - xl,){x - x^J, (22) 



with coefficients 



Al = m^, 



and discriminant 



Cl = T'ya-y)+mly + ml{l~y)-te, (23) 



Dl = BI-AAlCl. (24) 



Next, we introduce the following notation: = L*{x = y, y), L*|q = L*{x = 0, y) and 
^xy\y = ^xyi^ = y)- The two binomials are 



For Eq. fl20l) one has 



^*\y = "ni^ly + ml{l - y) - ie , 

kly\y = Q'y-ml{l-y). (25) 



dk'^ dV P'^^Tv ^Ti/ dL* 

dU = ^ = ^ ^ (26) 

{klyY ' dx LL* (L*)2 dx ' ^ ^ 



and the integral I{y) becomes: 

1 1 P'^dT 

m = HL'u) - ln(L| J] - - [ln(L%) - ln(Li J] jjy (27) 



After some more calculations we arrive at a one- dimensional integral over y, where the in- 
finitesimal parameter 6 cancels out. 
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2.4 Integration over y 

We proceed with the one-dimensional integrah 

1 

J= I dyl{y), 



with integrand 



and 



^2 



ln(L*|J-ln(p2^) 



2y/DT \y-yTi y- yr^ 



I,{y) = ln(C^)+ln(L*, )-21n(p2l/) 



Br 



Here P gets redefined: 



In 



'2CL + y{BL + VD^ 



Cr 



In 



'2CL + y{BL-y/Dl) 



Cr 



The quadratic trinomial Cl (see ( 123|) ). 

Cl = ary'^ + bry + ct, 

I. I 2 2 

Ot = J + - m4 , 
Cr = ml — ie , 



has the roots: 



bT±\/D- 



T 



2T2 



The following steps of integration with respect to y deviate from the presentation 
Ref. [3]. 

2.4.1 Splitting into three parts 

Let us redistribute terms in Eqs. (!29l) - (!30l) into three parts: 

I{y) = Io{y) + Iiiy) + l2{y), 



with 



h{yi 



1 1 

+ 



1 



2VDt yi 



Dt \y - z/Ti y- yT2 



In 1 - 



yio 



-In 



P^y — ie 



uia — le 



For hiyi) and /2(?/2) we have changed the variables, 

yi^y-yTi, 

and used the notation 



A^T^ + ml, 

-A^yn + A{T'' + A13) + 2m^ A34 , 
(r2-m^)yTi-T2-Ai3. 



Furthermore, 

C(yi) 

and 



'y2 



In 
In 



ml — ie 



y 



- In 1 - ^ 

Cjyi) - mljyi + yrj^ + m| - + (y^ + yrJ(A34 + y/Dl 

C{yi) 



yio 
yd 



A 



34 



= {Q"^ + ml){y - yd) , 



2 2 
mj — 



(35) 



(36) 



(37) 



(38) 



(39) 



These expression are still valid for all J's and are exact in m^. The main idea of this splitting 
arises from the observation that the 1714 singularities are completely confined to Iq, hence 
when calculating Ii and I2 one may take the limit 777.4 — before taking the integrals. 
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2.4.2 Part Jq 



The integral for Jq is straightforward, and we hmit ourselves to presenting the answer in 
limit 1714 — > 0: 



where 



+ 



^1 



-In 



P2 



In 



rp2 



-In 



p2' 

^ ) In {Rn) 



R 



T^ + mj- ie 

rr2 I 2 

i + 771]^ — 



Li2(l-i?i3) + C(2)L 



li 



, i = l,3 



and Li2{x) is the dilogarithm function defined by Li2{x) = — 



^ \n\l-xt\ 
t 



dt. 



2.4.3 Part Ji 

For Ji we directly take the limit m4 = in which it becomes 



Ji= h{y)dy = dy 



1 1 



In 



ms — le 



with the ratio 



where, since we neglect 777,4, 



Riv) 



T^3 



= D^{y) = Ay-2B,y + Ct, 

= A = T'^ + ml, 

= AT^.^-2T\ml-ie), 

= -{T^ + ml~ml), 

= + mj + ml - 2ie . 



Then we use the substitution: 



where 



E^C^^D^{y)\y=^-Cl + B^. 
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From Eq. fl45p we have: 



dx 



dy 



d{x) ' 
Cin{x) 
d{x) ' 
qEn{x) 
d^{x) 

dx 



Cix 



where 



d{x) 



n[x] 



A^Cf - (Ex - BiY , 



In order to take the integral (H2!) we need to know the hmits of variation of x. Let 

'5<1/<1, 

then 



A'Cf - Bf 
2CfE 



6 <x <1. 



Next, replacing the variable y in the ratio (143|) by x, we find 

{T^, + - Ex/{B^ + ACr)][l - Ex{Tf, - C^)/{Tf, + C{)/{B^ - AC^)] 



R{ 



XI 



{Tl, - C,)[l - Ex/iB, - AC,)][1 - Ex{Tl, + C,)/{Tl - C,)/{B, + AC,)] ' 



Taking account of all above equations, we rewrite the integral Ji of Eq. 
variables x oi y where convenient: 



using 



Ji 



1 



+ ln 



In 



+ - ie 



Tfg - Cl 



mg — le 



^ dx 

X 



'dy 

5 y 



'^ln(l 

y 



T'y 



+ mi — ie 



+ 



rp2 



dx 



X 



-In 1 



Ex 



+ ln 1 



Ex 
(fii - AC,) 



-In 1 



In 1 - 



EjTl - C,)x 
{Tf^ + C,){B,-AC,] 

E{Tf^ + Ci)x 
(Tl-C,){B, + AC,) 
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Finally we get 



T^ + mlJ '\Bi + ACiJ '\{Tl + Ci){Bi- ACi 

2.4.4 Part J2 

For 12(2/2) of ( l35l) we set 

y2 = y-yT2, (54) 

where in the limit — >■ 

2/T2 = 1 + ■ (55) 

• Transition from variable ?/2 to variable ^2 

At this step we make the transition from the variable ?/2 to ^2 by 



= ^jA^yl - 2^22/2 + C| = IC2I + 2/2^2 . (56) 

The presence of IC2I leads to the appearance of two branches, IC2I = ±6*2 > 0, in the final 
answer for J2. For our definitions see Eqs. (!38l) . 
As a consequence of (J56l) we get 

52 + ^2^2 ^ 



= 2 .2_.2 ' V ^2 = ^j— 2 , (57) 

J-^2 (-2 ^2 ''2 

with 

Ar2 = AIC2 + 2S2t2 + C2tl . (58) 
The Jacobian of the transition is 

^ - 2 (59) 



We also need 



dy2 ^ N^dh ^ ( C2 ^ ^ , 

7/2 (52 + C2t2)(^2_t2) + + ^ + ^^1 ^2, IDUJ 



and 

^7/2 (^^2 



2/2 -B2 + C'2t2 
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(61) 



The limits of integration are 



1 - yT2 



-yT2 



The remaining variables read 



max 



mm 



{m\ + Ai3)2 — Aml{ml — ie) , 
(T2 + Ai3)2 - Aml{ml - ie) . 



Replacement of variable from t2 to y' = y. 



The next replacement reads: 
where 

The Jacobian of this transition is 



t2 = A22/+tr, 



7-) jmax jmm 

J^t2 —^2 ~^2 ■ 



dy 



B2 + C2t2 C2{y -yd2)' 
The root yd2 is labelled by a second index ±, depending on the sign of C. 



y^2 



B2 + Ca^t^^" 



C^D 



2 ^t2 



For 12(1/2) we get from Eq. fl35l) 
1 



h{y2 



2^y2 



N{y2) + 



Do 



M{y2,VD2)-M{y2,-VD2) 



and 



= ±C2. 



The second term in curly brackets of fl68l) . after replacement of variables fl56|) . reads 



y2vD2 



1/2 



M(y2, V^)-M(|/2,-V^) 
M(t/2,/D^)-M(|/2,-v^) 



y2vD2 



C; 



1/2 



In 



In 



TT-n 



rir. 



In 



In 
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n, 



n. 



++/- 



-2T2(E2 + C2±t2), 

±2 (52 + C^t^ (t2 T A) , 
±2 {B<,^ACt) (t2T^). 



Solution C2 

h{y2) 



1 



with 



2VDTy2 



Co — Co 



N{y2) 



1 + 



Co 



In 



rir, 



do 



In 



2^2 



For and part of M we derive: 

_T2 



Niy2) 



In 
In 



n. 



-+ 



Ho 



d^ 

My 2) 



In 
In 
In 
In 



yi-Dt2 



In (1/2), 



+ In ( 1 - ^ 

Vi- 



2^2 

yi+Dt2{yi- - yd2) 
T'^yd2 

2C2?/d2 



+ In 1 - 



y 

yi+ 



+ ln 1 



y_ 

yd2 



-In 1 



- In ( 1 - 

yi- 



y_ 

yd2, 

In / 1 



y_ 

yi+ 



where 



yd2 
yi- 
yi+ 



B2 + t^^"C2 
Dt2C2 



Dt2 
4.m 
''2 



A + 1^*" 



D 



t2 



Putting the variable substitutions and Eqs. fl7^ into Eq. fl72|) . one gets: 



J2 



In 



ms — ze 



In 



1 - 1/T2 
-yT2 



+ J2, 



with 



1 



2 2(T2 + mf 

2C2yd2 



In 



+ 



In 



yi_yi+Dt2 



2C2yd2 
yi-yi+Dt2 
yi^Dt2 

J^2 



[nyi+) + iyi-)] + liHyi^) + liHyi^) 



-In 



+ ln 



yi+Dt2{yi- - yd2) 



d2 



I {yd2) 



+M{yi-,yi+) - Miyd2,yi+) + Miyi+,yi-) - M{yd2,yi-) - 2M{yi-,yd2) 
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and 

and the "master integral" : 
dy 



l{y) = In 1 



(7^ 



MiyuVd) 

• Solution C2 



In 1-1 

[y-Vd) V yi 



In 1 



yd 
yi 



Kvd) - Li2 



1 - z/d 
yi - yd 



Li2 



-yd 



yi - yd 



^2(2/2) 



1 



2VD^y2 



iV(^/2)^ 



with 



= C+ = - 1 + ^ 



In'"^ 



-In 



n6 



(79) 



^0) 



mr, + (mn — ie) . 



"2 — \ ^ ^ I ""2 ^ V"'3 

The derivation is similar and we limit ourselves to presenting the answer for 
2{T^ + mf)[ \yi-yi+Dt2j 2 2 



^1) 



+ 



In 



2C2yd2 



_ In I _y}+RR \ yi-Dt2{yd2 - yi+] 



I {yd2) 



yi-yi+Dt2j "' \ J \ T'^yd2 

+M{yi-,yi+) - M{yd2,yi+) + M{yi+,yi-) - M{yd2,yi-) -2M{yi+,yd2, 



12) 



2.5 Definitions of functions J 



thud 
sub 



The mass singularities in arguments of the logarithms may be compensated by combining J 
with one more Cq function: 



Coi-mi,-miQ';M^,md,0), (83) 



1 + 



ml + T2 



Co(-mi-m^,Q^;M^,m,,0). (84) 



where = U"^ for cross functions. The two mass-singular Cq functions appearing in Eqs. (183|) 
and (|84l) cancel in the total expression for the EW correction which proves the absence in it 
of logarithmic mass singularities (not KLN theorem!). 
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2.5.1 Calculation of the subtracted function J*!^^'^ 



The first function J*|^f is given by equation (|83|1 . where we may neglect m„ which does not 
give rise to a mass singularity: 



J: 



tbud 
sub 



J- 1 + 



r2 + m? 



Co(0,m4,(5 ;m3,m4,0) . 



For Cn we have: 



Co(0,m^,(5^;m3,m4,0) 



with 



In f-^) In(i?3) + Li2(i?3)-C(2) 



5) 



^6) 



^7) 



For Jf^ij^ one derives an expression for Jq that is free of the m4 mass singularity; Ji_2 remain 
unchanged: 

jthud jtbud I 7 I 7 fQQ\ 
•J sub - ^subfi + Jl + ■J2, lOOj 



jtbud 
'^subfl 



In - 



rrin 



-In (i?i3)-ln {Ru 



-Li2(i?3)-Li2(l-i?i3) L (89) 



2.6 Definitions of the function Jf^^^^i, 

If we want to neglect the h quark mass, m;,, we must perform a second subtraction of the 
mass singular Cq function Co{—mf, —ml, Q^, My^, rrib, 0) that appears in the limit mi, 0. 
Note that only one of J'^^'^ contains an mi, mass singularity. 

jfutubAQ^T';mi„mt,M^) = jf:,\Q\T^-mi„m,, M^) 



Q"^ + ml 
2^2 



Co{-ml-mlQ^-M^,mi„0). (90) 



Again, the mi, mass singular Cq function Co^—mf, —ml,Q'^; Myy,mii,0) cancels in the total 
EW correction. 

Since we do not want to consider the limit mt = 0, we simply rename the second function: 
Jl':rbU2iQ'^T';mt,m,,M„) = jfj'iQ' ,T';mt,m,, M^), (91) 
assuming mt = for this, non-singular case. 

2.6.1 The first function Jl'^'^ in the limit m^ ^ 

Here we simply present the limits of Jo,2 which develop mi mass singular logarithms; for Ji 
it is sufficient to set mi = in Eq. ( l53l) . 



J^2 



mn 



2 In 



T^-ie\ , /mi 



mo 



m 



- - In - 



mi 



In 



m^ 



- In ( ^ 1 In 

mi 



- Li2 (i?3) + Li2 

15 



-ml- ie 

2 

mi 



(92) 



Here P is given by Eq. fl31l) . For the J2 part one finds: 



J2 = ^ 



where 



and 



1 



Inf!^ 



+ 



21nl^l+ta 



y2 



mi 



In 



2 2 
m3m2 

^23 



+A^(y,_,y,+) - 2Li2 



T2 



A 



23 



A 



(r2-ze) (A23 + n 



yi+ 



m 



23 



2 ' 



^23 = ""^2 ~ ("^3 ~ "^e)- 

Summing Eqs. fl92p and fl93l) with Eq. flS^ in the hmit of mi = we get: 

'^iub'^ = Jo + Jl + J2- 



1 - y^- 



(93) 

(94) 
(95) 

(96) 



2.6.2 The function Jf-^,^, 

The Co function needed in Eq. (!90l) with exphcitly separated out m2 mass singular logarithm 
looks as follows: 



1 



+ 



In 



-Li2 



P2 



mi 



+ U2 



2m 



In 



m. 



m| 



In 



p2 



p\ + mf 



In 



(Q^ 



Pi + mf 
mf 



+ Lis 



P2 



+ mf 



2 2 
mfpf 



+ Li2 



-2m 

(P^)* 
Pi + mf 



In 



pI + mf 



2 

mf 



2C(2) . (97) 



Here P^ = + mf — ie and 



Let us emphasize that in this section Co has its own list of dummy arguments! 
In order to derive Jl^g^b 1 is sufficient to redefine Jq of ( |92l) into Jq by summing 



the first row of J2 of fl93|) and the Cq function of (1971) with the coefficient of Eq. 
collect together all rrib mass singular terms; that is: 

jtbud T' \ T \ 

'^subsub,! — ^0 "T ^1 "T "^2 5 

where Ji remains unchanged, given as before by the limit of Eq. fl53l) at mi 
denotes the rest of J2 without its first row. 



it with 
1, i.e. 



0, 



(98) 
and J2 



16 



For Jq we get: 



J'o 




with A23 defined by Eq. 

In Tabled] we give a comparison of real and imaginary parts of the function J^l^^f^b 1 defined 
by Eq. fl90|) and related ones, computed with the aid of the LoopTools package [5], vs numbers 
derived exactly from Eq. fl^JSl) with Jg given by Eq. ^% and Ji and by previous equations 
as explained just below Eq. (jHSD, i-e. results of this paper. 

The numbers are given for four values of s (the first two values near the kinematical edges 
and the latter two in the ±1 GeV vicinity of the W resonance) and for three values of cos^ 



at m. 



7 



10 



-40 



mi 



in GeV); first lines 



10 m2 = 174.3, ms 
LoopTools, second lines 



80.403, 7714 
- this paper. 



5-10 (all masses are given 



Table 1: Comparison of J\ 



cos 


= IGeV 


^/s = 173.2GeV 


-0.999 


-0.374523885975E-7, 0.483905777792E-7 
-0.374523675918E-7, 0.483905777797E-7 


0.220728294741E-1, -0.404470353912E-1 
0.220728294741E-1, -0.404470353912E-1 





-0.103142713149E-3, 0.143516783301E-3 
-0.103142713128E-3, 0.143516783301E-3 


0.490255707173E-1, -0.933810555325E-1 
0.490255707173E-1, -0.933810555325E-1 


0.999 


-0.779612000864, 1.57384208867 
-0.779612000864, 1.57384208867 


99.5783633460, -218.316883034 
99.5783633453, -218.316883034 




= 79.403GeV 


^ = 81.403GeV 


-0.999 


0.227686132083E-4, -0.303846732065E-4 
0.227686132074E-4, -0.303846732065E-4 


0.243884328473E-4, -0.327139389502E-4 
0.243884328481E-4, -0.327139389502E-4 





-0.504773282953E-3, 0.120287840705E-3 
-0.504773282955E-3, 0.120287840705E-3 


-0.507195719751E-3, -0.248876535390E-3 
-0.507195719749E-3, -0.248876535390E-3 


0.999 


-5.59686285604, 1.92566453330 
-5.59686285604, 1.92566453330 


-5.65305814019, -2.07873259553 
-5.65305814019, -2.07873259553 



As is seen from the Table, there is agreement within 9-12 digits for real parts and within 
12 digits for imaginary parts, which seems quite satisfactory given that we are using only 
Double Precision in the Fortran code. We also note that exact formulae with extremely small 
values of masses tti^, mi, m4 together with large masses m2, ma occur in LoopTools which 
could result in loss of computational precision. Unlike LoopTools, the formulae derived in this 
paper are rather compact and do not explicitly contain the masses m^, mi, m4; the latter 
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property is the main goal of this paper. Also, the short formulae of this paper demonstrate the 
underlying physics of singularity cancellations and, moreover, their execution is much faster 
than LoopTools. On the other hand, one may also conclude that this comparison proves the 
high reliability of the LoopTools package. 



2.6.3 The function J^^^Im 



This function is derived by means of simply setting m2 = in Eq. 
Since J'^^l) independent of m2, we have 



jtbud _ tO I tO I tO 
^subsub,2 — ^0 "T ^1 "T "^2 ' 



tO jtbud 

^0 ~ '^subfi 5 



with the latter given by Eq. 
For the part Ji one gets: 

J° = ln(/?i3)ln 



+Li2 



T2 + Ai3 + ie 



A 



r2 + 



13 



Lig 



A 



13 



R 



13 



+ 11, 



For the second part, J^, we introduce some common notation: 

1 



Let, 



A 



13 



(100) 



:ioi^ 



(102) 



2^2 



'Co 



— le. 



There are two solutions: 



2^2 



(103) 



if T"^ + ml-ml> 0, then 



1 



2(T2 + ml) 



where 



' 2^2 



yi^ 



In m^-- 1 - -In -y,^-- - In 1 - ^ 



2^2 



+l/^(,.,)-lnf^V(l/H) + Li2^ 



yi+Dt 



yi+ - yi- 



_j.2 _ 2(^1 - ie) 



Dt 
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yi- 



yu 



Lis 



-yi- 



yi+ - yi- 



i{yi 



(104) 



(105) 



• if T'^ + ml - mj < 0, then 



7° 




+ 2^ iyi+) + Li2 



(106) 



with 



Dt 

yi+ 



2T2 



(107) 



In Table [2] we give a similar comparison as in Table [U but now it is for the function Jf^ifg^i, 2 
defined by Eq. fllOOl) and related ones. The setup is the same as for Table [H Again, first lines 
— LoopTools, second lines — this paper. 



Table 2: Comparison of J\ 



cos 9 


^ = IGeV 


^ = 173.2GeV 


-0.999 


-0.943679024194E-3, 
-0.943679024191E-3, 


-0.495811769871E-6, 0.120460853637E-5 
-0.495811769847E-6, 0.120460853637E-5 





-0.177753225378E-3, 
-0.177753225357E-3, 


-0.373359924505E-6, -0.345101621985E-8 
-0.373359924465E-6, -0.345101621921E-8 


0.999 
SANC= 


-0.128202726721E-3, 
-0.128202726700E-3 


-0.253792663608E-6, -0.119084548438E-5 
-0.253792663600E-6, -0.119084548438E-5 




= 79.403GeV 


^/s = 81.403GeV 


-0.999 


0.917891464851E-3, 
0.917891464846E-3, 


0.843008928372E-3, 0.718400159039E-3 
0.843008928376E-3, 0.718400159039E-3 





-0.188317064868E-3, 
-0.188317064870E-3, 


-0.177932415951E-3, -0.901230697461E-4 
-0.177932415950E-3, -0.901230697461E-4 


0.999 


-0.266362647090E-3, 
-0.266362647091E-3, 


-0.258030469532E-3, -0.157505612576E-3 
-0.258030469531E-3, -0.157505612576E-3 



As seen from the Table, there is again agreement within 10-12 digits for real and imaginary 
parts. Below the W resonance this function is real. 
We switch now to another example of J functions. 
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3 Calculation of J for the process ud tb 

As in the previous section we shall consider only the direct and 'jW diagrams. For the process 
ud tb the list of arguments of the universal function J is: ml = mf, m\ = mf, m| = M^, 
ml = mq — singular mass. 

■J = ■JiwiQ'^^T'^'^'^d,mu,mb,mt) (108) 

In the following presentation we change the notation of masses as in ( 1131) and neglect the 
mass mu which does not lead to a mass singularity. 

3.1 Integration with respect to z,x 

We omit the details of the integrations with respect to z and to x. The standard ingredients 
are: 

N^y = -2k^yP2 = Ny , 

Ny = T\l-y) + Q^ + ml{l-y)+ml 

k^y — NyX — Ty , 

T^y = T^y{l-y)+ml{l-y) + mly, 
L = P^x-ie. (109) 

Contrary to the case of t — bud decay, the variables k^y and L* = L — kly are linear in x in 
this case since we have neglected m„: 

L* = L-kly = [-{T'' + ml){l-y) + ml-ml]x + T^ -ie. (110) 

3.2 Integration over y 

The one-dimensional integral with respect to ?/ — (1 — y) is 

1 

J= I dyl{y). (Ill) 







With the aid of Eq. (!27|) and with ingredients of Eqs. fll09H110"I) for the integrand I{y) one 
obtains a simpler result than Eqs. 



Ky\y Ty *e 



where 

)2 r^2 I 2 



Tl = T^y{l-y)+mly + ml{l-y) 
kly\y = mly{l-y) -mly + Q^{l-y) 

T*\y = -^Ivi^ - y) + ^ly + ^IC^ - y) - ■ (ns) 
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3.2.1 Integrand in the limit m4 — 

After some algebra one can derive the foUowing expression for I{y) valid in the limit —>■ 0: 



+ 



(T2 + ml) I 



'^{.y,yi) -'^{y,y2) 



Thus I{y) is expressed in terms of the auxiliary function 
1 



'^{y,v) 



(y - v) . 

with roots of the quadratic trinomial Ty — ie 

T'^ + ml - ie 

yi - 



In ( 1 - + In f 1 - - ln(l - - In 



y2 



The other objects arise from the trinomial k'^y\y in y 

-mly" + (A21 - g')2/ + (Q' + le), 



^y\y 
Bk 



A21 - 



Bl + Aml{Q^ + ie) 



with roots: 



2{-ml) 



yk2 



-Bk — y/Dk 



and trinomial L*^^ in y, i.e. 



A12 - mg, 



— AmKml — ie) 



with roots: 



yii 



2ml 



-Br, - .AD 



2ml 



3.3 Function J 

The final answer for the integral J of Eq. fillip is split into two parts, 

= ^0 + Jadd- 



(114) 



:ii5) 



(116) 



(117) 



:ii8) 



(119) 



(120) 



:i2r 
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The first part, free of mass singularity, is 



Jo 



1 



+Li2 



[ln(l -Z/fci) +ln(i?3)]/(?/fcJ -M{yLi,yki) -M{yL*,yki) 



1 -|/fei 



^M{yLi.yi) + Li2 } + {i/z.- - Z/lJ 



T2 + m? 

-[ln(l-|/i) + ln(i?3)]/(l/i)-Li2 



2/1 



1-2/1 



(122) 



and the additional part that depends on m\: 



J' 



add 



■ln(i?3)ln (i?i4). 



For simplicity we introduce the notation 

+ - ie 



(123) 



14 



(124) 



3.4 Subtracted function J^^f 

The subtracted function J^u^*** is defined by the following equation (see Ref. [1]): 



judtb T 



Co(0, Q ;m3,m4,0) . 



(125) 



After adding to this the mass singular Co function, we get the following expression free of 
mass singularity: 



J sub' = Jo + ln(-R3 



In - 



+ In {Rn) 



-Li2(i?3)+C(2), 



(126) 



where R3 and R13 are defined by Eqs. fl87|4ip respectively. 
3.4.1 Second substraction J^^bsubi 

A second substraction eliminates the h quark mass singularity: 



2^2 



Co{-mi,-mlQ';M^,mb,Q) . (127) 
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The limit in exists and reads 



udtb 
subsub,! 



1 



- In (Rs) In ( ^) + In 



P 



+2U2 



A 



23 



P 



with 



A23 

+ Li2 (1 + r23) - Li2 

rp2 



In 



:i28l 



A23 + ie 



- 3C(2) 



ri3 



mf 



TTin + ie 



(129) 



Tudtb 



See in this case Eq. f[T3|) for the meaning of rrii. 

In Table |3] we give the comparison of real and imaginary parts of the function J'^^bsub 1 
defined by Eq. (11271) and related ones, computed using the Loop Tools package with numbers, 
derived exactly from Eq. (11281) . i.e. with results of this paper. The numbers are given for two 
values of variable s (near threshold and at a high value), for three values of cos^. The masses 



are = 10 



-40 



nil = 10 



7712 



174.3, m3 = 80.403, m4 = 10 (all masses are given in 



GeV); first lines — LoopTools, second lines — this paper. 



Table 3: Comparison of J subsub 1 



cos 9 


^/s = 200GeV 


= lO^GeV 


-0.999 


-2.07194259988, 7.01180433443 
-2.07194259988, 7.01180433443 


-0.102517409116E-2, 0.102816191182E-2 
-0.102517409116E-2, 0.102816191182E-2 





0.878509046507E-3, 0.238029806504E-02 
0.878509046507E-3, 0.238029806504E-02 


0.895389809985E-6, 0. 184388800385E-6 
0.895389809985E-6, 0. 184388800385E-6 


0.999 


0.602158734688E-3, 0.913885117129E-03 
0.602158734687E-3, 0.913885117129E-03 


0.517213769582E-6, 0.486927343347E-7 
0.517213769582E-6, 0.486927343347E-7 



As is seen from this Table, there is agreement within 11-12 digits for the real parts and within 
12 digits for imaginary parts. We note again that Eq. ( I128P derived in this paper is very 
compact and does not contain the masses m^, mi, 7714 explicitly. 



3.4.2 No second subtraction Jsubsub2 

The function J subsub. 2iQ^ ^ f^b, Myy) with exchanged arguments has no mf, mass singula- 

rity. Since we need no limit in mj, we just rename the second function, assuming that m;, is 
set to zero: 



Jsubsub,2iQ', T'; 7nt, 7nt, M^) = J:TiQ', T'; m^, m,, M^) . 



(130) 



23 



The answer is found straightforwardly from Eq. 
interchanged arguments mi = m^, m2 = rrif, = 0: 



fr^r Tudtb 

I'J'^ ^sub 



judtb 



subsub,2 



1 



y2 

+Li2 
+Li2 



+ Li2 
-In 



13 



ml 



p2 
p2 



Li2 
hi 



A 



13 



p2 



m; 



Li 



and Jq from Eq. (11221) with 



mt 



A 



13 



It 



T'^ + ml 1 



m; 



l + ri3 



T2 + Ai3 

- Lia (Rs) - Li2 
- Li2 



T2 + m2 
T2 + mt 



In {R,)\n 



-C(2) 



T2 + m2 



mt 



ri3 



+ C(2) 



(131) 



The meaning of = remains unchanged. 

In Table H] we give a similar comparison as in Table [3] but now for the function J^^l^g^i^ 2 
defined by Eq. (11301) and the related ones. The setup is the same as for Table [H Again, first 
lines — Loop Tools, second lines — this paper. 



Table 4: Comparison of J\ 



cos 


^ = 200GeV 


^ = lO^GeV 


-0.999 


0.159231105282E-4, 0.118244938272E-3 
0.159231105282E-4, 0.118244938272E-3 


-0.195869068435E-3 0.240805881821E-3 
-0.195869068435E-3 0.240805881821E-3 





0.207059680614E-4, 0.881532449471E-4 
0.207059680614E-4, 0.881532449471E-4 


0.110155292282E-5 -0.167170803336E-6 
0.110155292282E-5 -0.167170803336E-6 


0.999 


0.234082370485E-4, 0.682484139707E-4 
0.234082370486E-4, 0.682484139707E-4 


0.459502689190E-6 -0.184916952415E-10 
0.459502689190E-6 -0.184916952411E-10 



As is seen from this Table, there is again agreement within 11-12 digits for real and imaginary 
parts. 



4 Functions J for the process hu td 

In this section we briefly consider J functions arising in four hu td processes. If one neglects 
the b quark mass, there appear only four doubly subtracted Jgulsub functions, see 

j's:lU2{Q^U\m„0,M„). (132) 

For four possible channels, one has the following correspondence between arguments Q^^T^, U"^ 
and s, cos 6, where always 6 = l{p2,P3) (see 4-momenta assignment below): 
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1. b{pi) + u{p2) C^(P3) + Kva) 

Q' = "-^c^ = -t, T'=-l^c^ = -u, U' = -s, (133) 

2. b{pi)+d{p2) ^u{Ps)+t{p4) 

Q2 = iZJ!ic_ = -t, T^ = -s, U^ = iZ^c+ = -u, (134) 



3. b{pi) + U{p2) -> rf(p3) + t(p4 



Q2 = ^-^c_ = -t, t2 = ^-!^c+ = -w, f/2 = -s, (135) 



4. 6(pi) + rf(p2) u{p^) + t(p4 



Q2 = = -t , T2 = -s f/2 = ^-^c+ = , (136) 



where 



c± = l±cos^. (137) 

There is an important difference between J functions considered in Sections [2] and [3] and 
in this Section. In the first two cases two Mandelstam variables = —t and f/^ = —u 
have the same sign; they are different only by the sign of cos^. In this, third case, has 
always the sense of — t, while and [/^ change their meaning —u or —s and therefore the 
sign. This is the reason why two t channel J's, for which both arguments are positive, can 
not be computed using the J functions derived for channels considered previously, and the 
calculation for such sign assignments has to be repeated from scratch. This task is beyond 
the scope of this paper and will be presented elsewhere. For the time being we will use a 
pragmatic solution. We noted that equations of Section [2] give correct answers for the case 
> and T2 < or < 0. 



Table 5: Comparison of J, 



cos 


= 200GeV 


^ = lO^GeV 


-0.999 


-0.173090211233E-04, 0.216758293310E-04 
-0.173090211232E-04, 0.216758293312E-04 


-0.467497992886E-06, 0.254446155152E-06 
-0.467497992836E-06, 0.254446155157E-06 





-0.451 130695567E-05, 0.149343040736E-05 
-0.451 130695569E-05, 0.149343040746E-05 


-0.462653095799E-06, 0.210944911427E-06 
-0.462653095748E-06, 0.210944911431E-06 


0.999 


0.172894787644E-04, -0.216698015838E-04 
0.172894787643E-04, -0.216698015840E-04 


0.105671318711E-06, -0.193311004466E-06 
0.105671318762E-06, -0.193311004461E-06 



This is illustrated by Tables [SHHl where we give a comparison of real and imaginary parts 
of the function Jgulf^^i and Jgubf~i,2 the t channel, Eq. (11321) . for the case > and 
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Table 6: Comparison of J, 



cos 9 


^ = 200GeV 


= lO^GeV 


-0.999 


0.108310963092E-03, 0.297113203560E-03 
0.108310963092E-03, 0.297113203560E-03 


-0.781547133959E-06, 0.303109233776E-06 
-0.781547133958E-06, 0.303109233776E-06 





0.288401492494E-04, 0.656706410374E-04 
0.288401492495E-04, 0.656706410373E-04 


-0.723035031957E-06, 0.259579787125E-06 
-0.723035031957E-06, 0.259579787125E-06 


0.999 


-0.803387429231E-04, -0.296822348620E-03 
-0.803387429230E-04, -0.296822348620E-03 


0. 100393594508E-06, -0. 166942822727E-06 
0. 100393594508E-06, -0. 166942822727E-06 



< or t/^ < using the same formulae and setup as for Tables [T] and [21 first lines — 
LoopTools, second lines — this paper. 

As is seen from these Tables, there is agreement within 10-11 digits for real and imaginary 
parts. Note that one can reach agreement to all visible digits with Real*16. 

As far as two t channel J's are concerned, for which both arguments are positive, we accept 
a temporary solution for the time being, noticing that previously computed J functions return 
correctly only the real parts. We recall that imaginary parts do not contribute at the one-loop 
level. The real parts are illustrated by Table [71 



Table 7: Comparison of J, 



cos 9 


yfs = lO^GeV 


^ = lO^GeV 




jbutd,+ 
subsub,! 


jbutd,+ 
subsub,2 


-0.999999 


-0.372108952046D+01 


-0.122105779653D-02 




-0.372108952003D+01 


-0.122105779653D-02 


-0.9990 


-0.122256223637D-02 


-0.256945712509D-03 




-0.122256223637D-02 


-0.256945712509D-03 





0.707331687150D-06 


0.122739516453D-05 




0.707331687129D-06 


0.122739516453D-05 


0.9999 


-0.155145904970D-06 


-0.149761802457D-06 




-0.155145905078D-06 


-0.149761802457D-06 


0.99999 


-0.515690518852D-06 


-0.496884568518D-06 




-0.515690519134D-06 


-0.496884568519D-06 



As is seen from these Tables, there is an agreement within 9-12 digits for real parts and 
that the agreement does not improve with Real*16 computations. 



5 SANC packages 

The numeric comparison with the LoopTools library presented in this paper can be verified 
with help of the SANC software packages. We have three packages related to the Jsubsub,i{2) 
functions for three channels. They are available to download from the web pages of SANC 
project [6]. 
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5ANC JAW packages. 

O 9/12/2009 SANC JAW butd vl.OO package (6.6 Kb tgz-fila) [last stable version] 
O 9/12/2009 SANC JAW tbud vl.OO package (5,2 Kb tgz-file) [last stable version] 
O 9/12/2009 SANC JAW udtb vl.OO package (5,1 Kb tgz-file) [last stable version] 
This packages are intended for calculation of JAW functions arising at the reduction of infrared divergent box diagrams. 

Figure 2: SANC packages download page 

Each package contains the following files: 

• main files for both Jsubsub,i{2) functions; 

• Jsubsub,i{2) source files; 

• utility functions; 

• Makefile; 

• README, INSTALL, LICENSE and other information files. 

6 Conclusions 

The invention and usage of a new class of functions J, relevant to the Passarino-Veltman 
reduction [Ij, has become a standard step in the chain of calculations in project SANC. 
Originally they were introduced in [3]. 

Obvious advantages and disadvantages of these functions are: 

1. final results are compact and evidently demonstrate the fundamentals of physics; 

2. the subtracted JwA,snh functions have no mass singularities; 

3. their compactness results in stable and very rapid calculations; 

4. the functions of this new class are linear combinations of the standard Passarino-Veltman 
Dq and Co functions and their explicit form depends on the concrete channel of a process; 
namely, they depend of the choice of channel because we have no universal expressions 
for the J functions, 

5. for the analytic calculations we have, nevertheless, a consistent method; 

6. it would be desirable to improve the way of the analytic calculations in order to obtain 
various channels by a simple rotation of their arguments, as is done in the Loop Tools 
package for Dq functions. 
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